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Chimera states are complex spatio-temporal patterns in which domains of synchronous and asynchronous 
dynamics coexist in coupled systems of oscillators. We examine how the character of the individual elements 
influences chimera states by studying networks of nonlocally coupled Van der Pol oscillators. Varying the 
bifurcation parameter of the Van der Pol system, we can interpolate between regular sinusoidal and strongly 
nonlinear relaxation oscillations, and demonstrate that more pronounced nonlinearity induces multi-chimera 
states with multiple incoherent domains. We show that the stability regimes for multi-chimera states and 
the mean phase velocity profiles of the oscillators change significantly as the nonlinearity becomes stronger. 
Furthermore, we reveal the influence of time delay on chimera patterns. 
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The investigation of coupled oscillatory systems is 
an important research field bridging between non¬ 
linear dynamics, network science, and statistical 
physics, with a variety of applications in physics, 
biology, and technology^ 2 l The analysis and nu¬ 
merical simulation of large networks with com¬ 
plex coupling schemes continues to open up new 
unexpected dynamical scenarios. Chimera states 
are an example for such intriguing phenomena; 
they exhibit a hybrid structure combining coex¬ 
isting domains of both coherent (synchronized) 
and incoherent (desynchronized) dynamics, and 
were first reported for the well-known model of 
phase oscillators®. In this paper, we investigate 
the influence of the local dynamics of the oscilla¬ 
tors upon the resulting chimera patterns. Using 
the Van der Pol oscillator, which is a model al¬ 
lowing for a continuous transition between sinu¬ 
soidal and strongly nonlinear relaxation oscilla¬ 
tions by tuning a single parameter, we show that 
multi-chimera patterns with multiple incoherent 
domains are promoted by increasing the nonlin¬ 
earity of the local oscillator dynamics. 


I. INTRODUCTION 

The last decade has seen an increasing interest in 
chimera states in dynamical networks 5 13 . It was shown 
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that they are not limited to phase oscillators, but can 
be found in a large variety of different systems including 
time-discrete mapd^l time-continuous chaotic models!^, 
neural systems-^S^, and Boolean network^!. Moreover, 
chimera states were found in systems with higher spa¬ 
tial dimension j 7 ! 9 ! 13 2 °KHI. Together with the initially 
reported chimera states, which consist of one coherent 
and one incoherent domain, new types of t hese pe cu- 
liar states having multiple incoherent regions 1 - 6 ’ 18 * 23 * ^, 
as well as amplitude-mediated- 6 27 j, and pure amplitude 
chimera and chimera death stated 28 were discovered. 

In many systems, the form of the coupling defines the 
possibility to obtain chimera states. The nonlocal cou¬ 
pling has generally been assumed to be a necessary con¬ 
dition for chimera states to evolve in coupled systems. 
However, rece nt stu dies have shown that even global all- 
to-all couplin;£ 2 -®2M]] as well as more complex coupling 
topologies allow for the existence of chimera states^" 37 . 
Furthermore, time-varying network structures can give 
rise to alternating chimera stated^. 

The important question of the main features that give 
rise to chimera states in coupled systems has been widely 
discussed, but no conclusive answer has been given yet. 
In systems of phase oscillators, the value of the phase lag 
parameter a , which occurs in the coupling function, is 
crucial. In nonlocally coupled systems, the range of the 
coupling and its strength play the key role. If the local 
dynamics of each unit is described by a two- or higher di¬ 
mensional system, then the interaction scheme between 
the units plays an important role, i.e., which variable is 
coupled to which variable of the other nodes. Chimera 
states have also been shown to be robust against inho¬ 
mogeneities of the local dynamics and coupling topol- 

0 ^) 3333 . 

Possible applications of chimera states in natural and 
technological systems include the phenomenon of uni- 
hemispheric sleepf®, bump states in neural system^® 2 !, 
power gridspA, or social systemd®! Many works consider- 
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ing chimera states have mostly been based on numerical 
results. A deeper bifurcation a nalysi s 45 and even a pos¬ 
sibility to control chimera stateP^ED were obtained only 
recently. 

The experimental verification of chimera s tates was 
first demonstrated in optical and chemica l 49 ! 50 * sys¬ 
tems. Further experiments invol ved mechanicaP^, elec- 
tronicP2l!31 and electrochemical^*^® oscillator systems, 
Boolean network^®!, the optical comb generated by a 
passively mode-locked quantum dot laser- 5 ®, and super¬ 
conducting quantum interference device^ 5 ®. 

In previous investigations of chimera states, usually 
the character of the local node dynamics has been con¬ 
sidered as fixed. In the current study, we address the 
issue of the impact of the local dynamics. We analyze 
the properties of chimera states, when the dynamics of 
individual oscillators smoothly changes from sinusoidal 
to nonlinear relaxation oscillations. For this reason, we 
choose the Van der Pol oscillator to describe the dynam¬ 
ics of each node. The Van der Pol oscillator^ has a long 
history of being used in both the physical and biological 
sciences, as a generic model for electrical circuits®*! and 
action potentials of neurons, respectively. 


II. THE MODEL 


In our study, we consider a system of nonlocally cou¬ 
pled Van der Pol oscillators with ring topology, where 
each element of the system interacts with a fixed range 
of its neighbors in both directions: 

u k =e(l - u\)u k - u k 
0 . j-k+R 

+ 2R Y - u k) + b 2 (u 0 - U k )}, 

j—k—R. 

with k = 1,..., N where all indices are taken modulo N , e 
is the bifurcation parameter of the individual oscillator, 
a denotes the strength of the coupling, R is the number 
of coupled neighbors (in each direction), and bi,b 2 are 
the interaction parameters. For such a form of coupling 
it is convenient to consider the ratio r = R/N , which we 
denote as a coupling range. The uncoupled Van der Pol 
oscillator has a stable trivial steady state u = 0 for e < 0 
and exhibits a supercritical Hopf bifurcation at e = 0. 
Here we consider e > 0. 

Introducing a new variable v k = u k , Eq. ([!]) can be 
rewritten in the form of a two-dimensional system: 


u k =v k 

v k =e(l - u 2 k )v k - u k 
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a 

2 R 


j=k+R 

Y [bi(uj 


j= k -R 


~u k ) + b 2 (vj 


v k)} ■ 


( 2 ) 


The form of the coupling in the system Eq. 0 or 0 is 
inspired from biological systems, describing interaction of 


the cells or pattern generation in locomotio d -— * — l A sim¬ 
ilar form of the coupling is also used in mechanics®^. The 
cross-couplings between the u- and the v- variable play an 
important role, they were shown to be necessary for the 
existence of chimera and multi-chimera states in systems 
of nonlocally coupled FitzHugh-Nagumo oscillators 16 . 

The dynamics of the system Eq. |2]) is determined by 
five parameters: e defines the dynamics of each individ¬ 
ual unit, and the parameters a , R , &i, and b 2 specifies 
the coupling. In order to find suitable values for some of 
the system parameters in the regime where Eq. 0 can 
describe chimera states, we will use the experience from 
simpler systems of coupled Kuramoto phase oscillators. 
For this reason we transform our system using the phase 
averaging technique on a rotating frame for slowly vary¬ 
ing amplitude r k and phase 8 k : u k {t) = r k (t) sin (t+9 k (t)) 
and v k (t) = r k (t)cos(t + 9 k (t)). As a result, we obtain 
the approximate system 

( 4 ~S (2jR+i)6a ) 

_ k+R 

+ JE\/ b i + b i Y r j cos{6 k -9 j +a) 

j—k—R (3) 

9 k ~(2R + l)b 1 

- k-\-R 

“^?V 6i + & 2 51 ^sin {Qk-Oj+a) 

j—k—R k 

with a = arctan(foi/& 2 ), b 2 > 0, and k = 1,..., N. 

The parameter a in the system © can be associated 
with the phase lag parameter in the systems of coupled 
phase oscillator^. This parameter is crucial for the ap¬ 
pearance of chimera states in the phase oscillator net¬ 
work. IiP it was shown that a value of the phase lag 
parameter close to but slightly less than tt/2 allows for 
the existence of chimera states. 

In the following, using the experience from the phase 
oscillator network, we fix the interaction parameters to 
be b\ = 1 and b 2 = 0.1, such that a « 1.47 is close to 
7t/2. With this parameter choice, we will focus further 
on the original system Eq. ([2]) , and vary the parameter e 
that defines the type of local dynamics of each element, 
as well as the coupling parameters a and r describing the 
strength and the range of the coupling, respectively. 



III. THE IMPACT OF LOCAL DYNAMICS 

Varying the bifurcation parameter e results in a change 
of the character of the local node dynamics. If e is small, 
the uncoupled individual elements of the system perform 
harmonic oscillations on a limit cycle, which is approx¬ 
imately a circle. With increasing e, the individual limit 
cycle becomes distorted and changes its form to relax¬ 
ation oscillations. 

Figure [T] demonstrates examples of chimera states for 
the system of N = 1000 elements, e = 0.2, and decreas- 
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FIG. 1. Snapshots of the variables Uk (upper panels), mean phase velocities oJk (middle panels), and snapshots in the phase 
space (uk,Vk) (bottom panels, limit cycle of the uncoupled unit shown black), (a) r = 0.35, a = 0.05, (b) r = 0.2, a = 0.09, 
(c) r = 0.13, a = 0.09, (d) r = 0.1, a = 0.09. Other parameters: N = 1000, 6 i » 1, 62 = 0.1, and e = 0.2. 




FIG. 2. Stability regimes for multiple chimera states, (a) e = 0.2, black squares marked by A-D show parameter values 
corresponding to panels (a)-(d) in Fig. [I] The insets show snapshots of coherent spatial profiles for parameter values A*(r = 
0.35, cr = 0.12; I\ = 0), B*{r = 0.2, a = 0.28; A' = 1), C*(r = 0.13, a = 0.34; A' = 2), D*(r = 0.1, a = 0.38; A' = 3); (b) e = 0.4. 
Other parameters as in Fig. [l] 


ing coupling range. The upper panels depict snapshots 
of the variables Uk for fixed time T = 50000. As ini¬ 
tial conditions we use randomly distributed phases on 
the circle u 2 + v 2 = 4, i.e., around the limit cycle of the 
uncoupled system, which is approximately a circle of ra¬ 
dius 2. One can clearly distinguish coherent and incoher¬ 
ent domains, a characteristic signature of chimera states. 
Elements that belong to the incoherent domain are scat¬ 
tered along the limit cycle, as shown with red points in 
the bottom panels of Fig. [l] where the black line denotes 
the limit cycle of the uncoupled unit with corresponding 
value £ = 0.2. The individual nodes perform a nonuni¬ 


form rotational motion, but neighboring oscillators are 
not phase-locked. To illustrate this, the middle panels of 
Fig- 0 show the mean phase velocities for each oscillator 
calculated as Uk = 2nMk/ AT, k = 1 ,..., N, where Mk is 
the number of complete rotations around the origin per¬ 
formed by the fc-th node during the time interval AT. 
Throughout the paper we use AT = 50000 for the calcu¬ 
lation of the mean phase velocities cok corresponding to 
several thousand periods. The values of uik lie on a con¬ 
tinuous curve and the interval of constant uik corresponds 
to the coherent domain, where neighboring elements are 
phase-locked. This mean phase velocity profile is a clear 
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indication of chimera states and similar to the case of 
coupled Kuramoto phase oscillators 3 4 . 

In addition to chimera states with one incoherent do¬ 
main [Fig. |TJa)] , we observe chimera states with multiple 
incoherent domains shown in Fig. 0b) -(d), i.e., multi¬ 
chimera states. In the following, we will use the notation 
n-chimera for a chimera state with n coherent and n inco¬ 
herent domains. In analogy with networks of phase oscil¬ 
lators, here we observe chimera states where the number 
of incoherent domains is even. The number of incoherent 
domains increases with decreasing coupling range. 

Figure [2|a) shows the stability regimes for chimera 
states with one and multiple incoherent domains in the 
plane of coupling range r and coupling strength ct for 
e = 0.2. Indeed, for large coupling range, we observe the 
stability regime for chimera states with one incoherent 
domain, and regimes for chimera states with two, four, 
and six incoherent domains follow subsequently with de¬ 
creasing coupling range. The overlaps of these regimes 
are characterized by multistability, when each of the 
chimera states can be obtained in the system depend¬ 
ing on the choice of the initial conditions. The regimes 
shown in the diagram are obtained by starting from the 
chimera states shown in Fig. [l] and using this pattern as 
initial condition for the neighboring parameter set, and 
so forth with a step size of Ar = 0.01 and Act = 0.01. 
Black squares denoted by A-D show values of the param¬ 
eters (r, ct) that correspond to the examples presented in 
Fig. [lja)-(d), respectively. 

For larger coupling strength ct, we observe coherent 
states in the system ([2|. They are characterized by the 
wavenumber K defining the number of maxima (min¬ 
ima) in the spatial profile, and K = 0 corresponds to 
complete in-phase synchronization. The wavenumber in¬ 
creases with decreasing coupling radius, and exemplary 
snapshots are shown in the insets of Fig. [2](a). For large 
coupling strength ct, the system is characterized by high 
multistability, and depending upon initial conditions one 
can obtain coherent solutions with different wavenum¬ 
ber. In our system there exist two different types of 
chimeras, amplitude-mediated chimeras and pure phase 
chimeras, and these are generated by different bifurca¬ 
tion mechanisms. The amplitude-mediated 2-, 4-, 6- 
chimeras (Fig. 2(a), snapshots in Fig. l(b),(c),(d), top 
panel) are generated from smooth, completely coherent 
spatial profiles of wavenumbers K = 1,2,3, respectively, 
by a coherence-incoherence bifurcation with decreasing 
coupling strength as indicated by the insets in Fig. 2(a). 
At the onset of chimeras the smooth coherent profiles 
break up into spatially coherent domains corresponding 
to the upper and the lower parts of these profiles, and 
incoherent domains in between. Therefore these incoher¬ 
ent domains occur in pairs (2-, 4-, 6-chimeras). Such 
coherence-incoherence bifurcations have also been ob¬ 
served for other local dynamics, ag., logistic maps and 
Rossler system^, the cosine mapP^, and Stuart-Landau 
oscillator^!!. In contrast, the pure phase chimeras (1- 
chimera in Fig. 2(a),(b)) arise from completely in-phase 


synchronized (K = 0) profiles. The shift of the regimes 
for (multiple) n-chimeras to smaller coupling range r 
with increasing n is typical for various nonlocally cou¬ 
pled systemsPH!£128l48J 

For small e, the limit cycle of each individual Van der 
Pol oscillator is close to a circle, corresponding to si¬ 
nusoidal oscillations, and the similarities to the chimera 
states in a system of phase oscillators are clearly revealed. 
However, the hybrid solutions we observe in the system of 
coupled Van der Pol oscillators demonstrate chimera be¬ 
havior both for phases and amplitudes. This can be seen 
in the bottom panels of Fig. |T] where red dots denot¬ 
ing the snapshot of all nodes deviate in their amplitudes 
slightly from the limit cycle of the uncoupled unit (black 
line). 

Figure [2](b) depicts the stability regimes for chimera 
states in the system ([2]) with e = 0.4. Compared to the 
case of e = 0.2, the regimes for the chimera states with 
multiple incoherent domains become larger, and chimera 
states in this case can be obtained for a wider range 
of coupling strength ct. The reason why the regions in 
Fig. |(b) (e = 0.4) are larger than those in Fig. [^a) 
(e = 0.2) is related to the following qualitative argu¬ 
ment: If the coefficient £ of the nonlinear term in Eq. (1) 
is increased, the coefficient ct of the coupling term has to 
be scaled up accordingly to balance the nonlinear term. 
Hence as e is increased the chimera regions extend to 
larger values of ct. 

When the parameter e is increased, the limit cycle of 
the individual uncoupled Van der Pol oscillators deforms, 
and the dynamics on the cycle becomes of slow-fast type. 
Further increasing £ leads to strongly nonlinear relax¬ 
ation oscillations. This will be discussed in the following. 

Figure [3](a) depicts the stability regimes for system K 
with £ = 0.8. Compared to the cases of smaller e, there 
are several qualitative differences in the stability regimes 
of chimera states. First, chimera states can be observed 
for a much larger range of coupling strength ct. Second, 
chimera states with one incoherent part cannot be ob¬ 
served in the system any more: for large coupling range 
we observe chimera states with four incoherent parts, and 
furthermore with decreasing coupling range the multi¬ 
plicity of the incoherent domains of the chimera states 
increases. Black triangles denoted by E, F, and G show 
values for the parameter pairs ( r , ct) that correspond to 
examples of chimera states depicted in Fig. Ka)-(c), re¬ 
spectively. 

The peculiarity of the diagram presented in Fig.^a) is 
the presence of two separate regimes for chimera states 
with four incoherent domains. Analyzing this diagram 
in more detail, one can see that there are two quali¬ 
tatively different regions. The first region appears for 
large coupling strengths, and contains stability regimes 
for chimera states with two and four incoherent domains 
(yellow region containing point E and blue region). These 
states can be characterized by strong amplitude dynam¬ 
ics, and the maximum values of the mean phase velocity 
profile correspond to the coherent domains of chimera 




5 


2-chimera 

4-chimera 

6-chimera 

8-chimera 




FIG. 3. Stability regimes for multiple chimera states, (a) e = 0.8, black triangles marked by E-G show parameter values 
corresponding to the panels (a)-(c) in Fig. [d] (b) e = 1.5, black circles marked by H-J denote parameter values corresponding 
to the panels (a)-(c) in Fig. [5] Other parameters as in Fig. [l] 








FIG. 4. Snapshots of the variables Uk (upper panels), mean phase velocities uJk (middle panels), and snapshots in the phase 
space ( Uk,Vk) (bottom panels, limit cycle of the uncoupled unit shown black), (a) r = 0.17, a = 0.8, (b) r = 0.35, a = 0.3, 
(c) r = 0.25, a = 0.2. Other parameters as in Fig. ga). 


states. The example presented in Fig. ga) [correspond¬ 
ing to point E in Fig. ga)] depicts these features. Com¬ 
pared with Fig. [l] for small e, we notice that the chimera 
states shown there also show distinct variations along 
the limit cycle, and the coherent domains of the chimera 
states correspond to the maximum in the mean phase 
velocity profiles. 

The second, qualitatively different part of the diagram 
in Fig. ga) includes three regions for small coupling 
strengths (yellow including point F, green, and gray). 
These regions form a similar sequence with increasing 
multiplicity of the chimera starting from four incoher¬ 
ent parts. However, they exhibit a qualitative difference. 


Inspecting Fig. gb) and (c), which show examples that 
correspond to the parameter pairs (r, a) denoted by F 
and G, one can notice that the amplitude dynamics be¬ 
comes weaker in these cases, and the network solution is 
close to the limit cycle of the uncoupled node shown as 
black line in the bottom panels. Moreover, the minimum 
of the mean phase velocities profiles now corresponds to 
the coherent domains of the chimera states. 

The difference between Fig. B«o on one hand, and 
Fig. gb) and (c) on the other hand is due to two 
different types of chimeras. The 2-, 4-, 6-chimeras 
in Fig. ga),(b) and Fig. ga) (point E), and the 
chimeras in Fig. gF,G,H,I,J) belong to two differ- 
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FIG. 5. Snapshots of the variables Uk (upper panels), mean phase velocities u>k (middle panels), and snapshots in the phase 
space ( Uk,Vk ) (bottom panels, limit cycle of the uncoupled unit shown black), (a) r = 0.4, a = 0.1, (b) r = 0.22, a = 0.1, 
(c) r = 0.17, a = 0.1. Other parameters as in Fig. 0b). 


ent types of chimeras: Fig. 0a), (b) and Fig. 0a), 
point E (corresponding to phase portraits shown in 
Fig. 0b), (c) ,(d) and Fig. 0a)) correspond to amplitude- 
mediated chimeras with strong amplitude-phase cou¬ 
pling, whereas Fig. 0a), points F,G and Fig. 0b) (cor¬ 
responding to phase portraits shown in Fig. 0b), (c) and 
Fig.0a),(b),(c)) correspond to pure phase chimeras sim¬ 
ilar to the ones found for Kuramoto phase oscillators, 
and since the phase oscillator model can generally be ob¬ 
tained from amplitude-phase models in the weak coupling 
limit, they occur in the stability diagram (Fig. 0a)) only 
for small coupling strength (points F,G), as opposed to 
the amplitude-mediated chimeras (point E). This differ¬ 
ence is visible in the phase portraits of Figs. 1,4,5 (bot¬ 
tom panels), where the spread of the various oscillators 
around the cycle of the uncoupled oscillator (black cycle) 
is large for amplitude-mediated chimeras, and very small 
for pure phase chimeras where the phase of the cycle is 
the only dynamical degree of freedom. The difference 
also shows up in the smaller amplitude variation of the 
mean phase velocity in the middle panels of Fig. 0b), (c) 
(pure phase chimeras) as compared to Fig. B a ) , and in 
the inverted tOk profiles: the coherent regions correspond 
to the minima (Fig. 0b), (c)) and maxima (Fig. 0a)), 
respectively. 

Further increase of the parameter e of individual Van 
der Pol oscillators leads to an even stronger deformation 
of the limit cycle. In the (r, a) parameter plane the sta¬ 
bility regimes for chimera states with four and more in¬ 
coherent domains can be observed as shown in Fig. 0b) 


for e = 1.5. The effect of the coexistence of two qual¬ 
itatively different types of chimera states is not present 
there any more, in contrast to the case of e = 0.8. Only 
the second type of the chimera states is observed in the 
systems now, and the stability regimes are enlarged to¬ 
wards larger coupling strengths. 


Figure 0 depicts examples of multi-chimera states that 
correspond to the parameter pairs (r, a) denoted by H,I, 
and J in Fig. 0b). The coherent domains of the chimera 
states correspond to the minimum of the mean phase ve¬ 
locity profile, and all oscillators stay very close to the 
limit cycle of the single uncoupled unit, thus the am¬ 
plitude dynamics of the chimera states in the systems 
with large £ is not as pronounced as in the networks with 
small e. 


We conclude that the nonlinearity of the local dy¬ 
namics indeed strongly influences chimera states in sys¬ 
tem (0. The character of the amplitude dynamics, the 
frequencies of the oscillators belonging to the coherent 
and incoherent domains of the chimera states, i.e., the 
mean phase velocity profiles, and the stability regimes 
in the coupling parameter plane undergo a qualitative 
change with variation of the parameter e. Stronger non¬ 
linearity (larger e) results in the dominance of multi¬ 
chimera states with weak amplitude dynamics. 
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FIG. 6 . Space-time plots of Uk for different values of time delay: (a) r = 1, (b) r = 3, (c) r = 6 . Other parameters: N = 1000, 
bi = 1, 62 = 0.1, r = 0.4, a = 0.1, e = 1.5. Initial conditions as shown in Fig. |5|a). Transients of 2000 time units are skipped. 
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FIG. 7. Snapshots of the variables Uk (upper panels), mean phase velocities tOk (middle panels), and snapshots in the phase 
space (wfc, Vk) (bottom panels, limit cycle of the uncoupled unit shown black), (a) r = 1, (b) r = 3, (c) t = 6 . Other parameters: 
N = 1000, bi = 1, fe 2 = 0.1, r = 0.4, a = 0.1, e = 1.5. 


IV. TIME-DELAYED COUPLING 

Together with the character of the local dynamics, the 
coupling between the individual units plays an impor¬ 
tant role for the properties of the chimera states. Time- 
delayed coupling if compared to the instantaneous one 
represents a more realistic way to model the interaction 
between the coupled units. Usually, the coupling range 
and strength influence the multiplicity of coherent do¬ 
mains in chimera states. However, it has been shown for 
phase oscillator networks that time delay can also induce 
multi-chimera^? 1 . The existence of chimera states in sys¬ 


t ems w ith time-delayed couplings has been also reported 
i] j 50 | 63 H 66] j n particular, for coupled phase oscillator sys¬ 
tems it has been found that chimeras are robust to small 
time delays and delay distribution^ and can become 
unstable depending on the value of delay®!. Here we 
consider a model that includes not only phase but also 
amplitude dynamics, and show how time delay in the cou¬ 
pling affects chimera states that exist in the undelayed 
system. We demonstrate that by varying the delay value 
one can both conserve and eliminate chimera patterns. 

Let us consider Eq. ([ 2 ]) modified by time-delayed cou- 





































































pling: 

Uk(t) =v k (t) 

v k (t) =e[l - u 2 k {t)\v k {t) - u k {t) 
j= k +R 

+ { & i [uj(t~ T ) ~u k (t)} + 

j=k—R 

+b 2 [vj(t-r) -v k (t)]} 

wth k = 1 ,N modulo N, where r is the delay time. 

Using the chimera state with four incoherent domains, 
shown in Fig.[5ja), as initial condition, i.e., as the history 
in the interval [—r, 0], we fix all system parameters cor¬ 
responding to this solution, and show exemplary space- 
time patterns of Eq. Q for different time delays. The 
period of a single uncoupled oscillator is close to 2n, and 
we neglect the transients of 2000 time units. For small 
time delay (r = 1), we observe a coherent traveling wave 
solution shown in Fig.^a). When the time delay is close 
to half the oscillation period (r = 3), the chimera pat¬ 
tern is stable and we continue to observe a chimera state 
with four incoherent domains as shown in Fig. |6jb). A 
larger time delay (r = 6), which is close to the period of 
a single oscillator, leads to complete synchronization of 
all oscillators, see Fig. [6^c). Figure [7] depicts snapshots 
of the variable u k (upper panels), the same snapshots 
in the phase plane ( u k , v k ) together with the limit cy¬ 
cle of uncoupled oscillator, and the corresponding mean 
phase velocity profiles (middle panels), for the solutions 
shown in Fig. [0} The explanation for the effect of delay is 
that the delay time interacts with the intrinsic timescale 
(oscillation period) giving rise to resonance phenomena 
as found generally for delayed feedback control of steady 
states, deterministic limit cycles, and noise-induced oscil¬ 
lations, if the delay is an integer multiple or a half-integer 
multiple of the intrinsic timescal^3. Delay has a favor¬ 
able effect on chimeras if r is a half-integer multiple, and 
a favorable (stabilizing) effect on the synchronized oscil¬ 
lations if it is an integer multiple, and may induce trav¬ 
eling waves if it fits with neither condition. In case of 
the chimera (Fig. [7](b)), the delay leads to much longer 
transients, so that with the same length of time inter¬ 
val used for the calculation of the mean phase velocity 
as without delay, the profiles are more smeared out, but 
qualitatively similar. 

Our numerical evidence shows similar results for other 
values of the bifurcation parameter e. For chimera states 
with one or two incoherent domains and sinusoidal char¬ 
acter of the oscillations, small delay can lead not only to 
traveling wave solutions, but also to chimera states with 
higher number of incoherent domains. 

These examples demonstrate that time delay intro¬ 
duced in the coupling can either suppress or preserve the 
chimera patterns depending upon the value of the delay 
time relative to the intrinsic oscillation period. 


V. CONCLUSION 

In the current study we have demonstrated how the 
character of the local oscillator dynamics influences 
chimera states in networks of nonlocally coupled Van der 
Pol oscillators. Changing the bifurcation parameter of 
the single oscillators allows us to interpolate continuously 
between sinusoidal and strongly nonlinear relaxation os¬ 
cillations. We have shown that nonlinearity facilitates 
multi-chimera states. 

For small values of the bifurcation parameter e (sinu¬ 
soidal oscillations) chimera states are characterized by 
lower multiplicity and more pronounced amplitude dy¬ 
namics, and the maxima in the mean phase velocity 
profiles correspond to the coherent domains. Moving 
towards the relaxation oscillation regime, with increas¬ 
ing e, leads to a higher multiplicity of chimera states, 
but weaker amplitude dynamics. In contrast to the previ¬ 
ous case, the coherent domains correspond to the minima 
of the mean phase velocity profiles, i.e., the profiles are 
flipped. We have also found that time delay in the cou¬ 
pling strongly affects the chimera patterns in the system, 
and can lead to chimera suppression and the formation 
of traveling waves and complete synchronization. 

We have presented (multi-) chimera states of different 
type: (i) pure phase chimeras, which are similar to those 
found for Kuramoto phase oscillators or weakly coupled 
amplitude-phase models, and (ii) amplitude-mediated 
chimeras with strong amplitude-phase coupling. 

Our findings give new insight into the intriguing phe¬ 
nomena of chimera states, and demonstrate that the 
character of the local dynamics has a strong influence 
on the chimera patterns in the whole network. These 
results could be useful from the point of view of applica¬ 
tions dealing with different kinds of oscillators, as they 
can be realized, e.g., in electronic circuits. 
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